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ABSTRACT 


Using  a  number  of  Geiger-Mueller  tubes  in  the  University 
of  Iowa  instruments  on  Pioneer  10  and  Pioneer  11,  we  have  C  ,  ~ 
measured  the  interplanetary  intensity  of  cosmic  rays  j(E^  >  80  MeV) 
over  a  complete  solar  activity  cycle  and  out  to  a  heliocentric 
distance  of  over  32  AU.  Free  space  values  of  the  background 
corrections  due  to  the  radioisotope  thermal  generators  were 
determined  during  passage  of  Pioneer  11  under  the  dense  rings  of 
Saturn  during  the  1  September  1979  encounter.  These  improved 
corrections  and  ones  inferred  therefrom  for  Pioneer  10  have 
supplanted  earlier  estimates  and  have  been  applied  to  all  data. 

During  the  period  1972-198*+,  intensity  vs  time  curves 
have  a  flat,  more-or-less  constant  maximum  during  1975,  1976,  and 
1977;  a  brief  (~  8  months)  minimum  in  mid-1981;  and  another  brief 
minimum  centered  around  early  January  1983-  The  ratio  of  maximum 
to  minimum  intensity  (normalized  to  same  radial  distance)  is  «  2.7. 
The  temporal  correlation  of  intensity  at  P10  with  that  at  Pll  is 
good,  after  allowance  for  propagation  delay,  as  is  that  with  data 
from  the  Deep  River  Neutron  Monitor.  The  max/min  intensity 
ratio  for  the  latter  is  «  1.25-  As  noted  by  us  and  others  in 
previous  studies,  the  decline  in  intensity  after  the  prolonged 
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maximum  occurred  primarily  by  an  irregular  succession  of  step¬ 
wise  decreases  with  subsequent  slow  recovery  such  that  successive 
decreases  overlapped  cumulatively.  The  logarithms  of  the  ratios 
of  the  counting  rates  of  a  GM  tube  on  P10  to  that  of  a  matched 
GM  tube  on  Pll  (after  allowance  for  propagation  delay)  were 
plotted  against  the  differences  in  heliocentric  radial  distance 
Ar.  Three  such  independent  plots  were  made,  one  for  each  of  a 
matched  pair  of  tubes.  The  three  least-squares-fitted  trend 
lines  through  the  entire  body  of  data  for  each  pair  of  detectors 
yield  integral  radial  gradients  G  of  +2.00  (±  0.07),  +2.07  (±  0.05), 
and  +2.12  (±  0.0 6)  percent  per  AU.  Our  adopted  mean  value  with 
estimated  uncertainty  for  all  known  causes  is 

G  =  +2.06  (±  0.20)  percent  per  AU 

for  the  radial  range  1  <  r  <  32  AU.  The  reciprocal  of  G,  namely 
50  AU,  is  the  scale  length  of  the  heliosphere.  By  direct  observa¬ 
tion,  the  modulating  influence  of  the  sun  is  still  clearly 
present  at  32  AU.  Hence  the  radius  of  the  heliospheric  modula¬ 
tion  region  r  >  32  AU.  Various  arguments,  no  one  of  which  is 

D 

convincing,  suggest  that  r  lies  between  50  and  100  AU. 

D 
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I.  INTRODUCTION 

This  paper  gives  a  progress  report  on  our  long  term  study 
of  the  interplanetary  intensity  of  the  cosmic  radiation  out  to 
great  distances  from  the  sun  using  a  system  of  detectors  on  the 
spacecraft  Pioneer  10  and  Pioneer  11.  The  observations  give  the 
integral  intensity  of  cosmic  rays  having  energies  sufficient  to 
penetrate  the  physical  shields  of  the  detectors,  -which  correspond 
to  the  range  of  a  proton  of  kinetic  energy  Ep  =  80  MeV  or  of  a 
helium  nucleus  of  specific  kinetic  energy  ea  =  80  MeV/nucleon.  The 
omnidirectional  integral  intensity  is  proportional  to  the  count¬ 
ing  rate  of  a  detector. 

As  of  1  April  1984,  the  heliocentric  distances  of  the  two 
spacecraft  were  52.48  and  15.81  AU,  respectively.  The  Jupiter 
encounter/heliospheric  mission  of  Pioneer  10  began  with  launch 
of  the  spacecraft  on  5  March  1972.  Following  its  encounter  with 
Jupiter  in  early  December  1975,  the  spacecraft  has  been  on  a 
hyperbolic,  escape  trajectory  from  the  solar  system.  The 
asymptote  of  this  trajectory  is  at  =  85?5,  3^  =  +2?9 
(heliocentric  ecliptic  longitude  and  latitude  in  1950.0 
coordinates)  and  v^  =  2.59  AU  yr  The  projected  heliocentric 
distances  r  during  the  next  ten  years  are  as  follows: 
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35  AU,  March  1985;  40  AU,  December  1986;  45  AU,  November  1988; 

50  AU,  September  1990  ;  55  AU,  August  1992;  and  60  AU,  July  1994. 

The  Jupiter  encounter/Saturn  encounter/heliospheric 
mission  of  Pioneer  11  began  with  launch  on  6  April  1973.  This 
spacecraft  encountered  Jupiter  in  early  December  1974  and  Saturn 
in  late  August  -  early  September  1979-  Since  the  latter  en¬ 
counter  Pioneer  11  has  also  been  on  a  hyperbolic,  escape  tra¬ 
jectory  from  the  solar  system  with  =  290.°2,  =  +12?8,  and 

v^  =  2.20  AU  yr”\  Projected  heliocentric  distances  are  as 
follows:  20  AU,  January  1986;  25  AU,  January  1988;  30  AU, 

January  1990;  35  AU,  January  1992;  and  40  AU,  February  1994. 

The  University  of  Iowa  instruments  on  both  spacecraft  con¬ 
tinue  to  operate  flawlessly.  As  Pioneer  11  flew  under  the  rings 
of  Saturn  on  1  September  1979  there  was  a  unique  opportunity  to 
determine  the  free-space  background  counting  rates  of  the 
detectors  caused  by  the  Radioisotope  Thermal  Generators  (RTG’s) 
that  furnish  electrical  power  to  the  spacecraft.  The  present 
paper  is  the  first  of  our  series  of  papers  on  the  interplanetary 
cosmic  ray  intensity  [Van  Allen,  1972;  Thomsen  and  Van  Allen, 
1974;  Van  Allen,  1976,  1979>  1980]  to  use  the  revised  background 
rates  derived  from  observations  during  the  Saturn  encounter. 
Hence,  this  paper  provides  an  improvement  over  previous  ones  as 
well  as  an  extension  of  the  study  to  greater  heliocentric-  dis¬ 
tances.  (See  Appendix  I  for  Remarks  on  Validity  and  Accuracy 


of  Data .  ) 


Relevant  recent  papers  using  data  from  various  nstruments 
on  Pioneers  10  and  11;  on  the  two  other  outer  heliospheric 
spacecraft  Voyagers  1  and  2;  and  on  IMP  7,  IMP  8,  and  Helios  1 
and  2  in  earth  orbit  and  in  the  inner  heliosphere;  and  from 
earth-based  neutron  monitors  are  by:  Webber  and  Lockwood  [1981], 
McDonald  et  al.  [1981a];  McDonald  et  al.  [1981b],  McKibben  et  al. 
[1982],  Lockwood  and  Webber  [198*0,  Venkatesan  et  al.  [1964a], 
Decker  et  al.  [1984],  and  Venkatesan  et  al.  [1984b]. 

The  trajectories  of  Pioneer  10,  Pioneer  11,  Voyager  1, 
and  Voyager  2  are  shown  in  Figure  1. 
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II.  OBSERVED  DATA 


Daily  mean  counting  rates  of  our  three  matched  pairs  of 
detectors  are  shown  as  a  function  of  time  in  Figures  2,  3,  and  4 
Occasional  upward  spikes  are  caused  by  solar  energetic  particles 
All  counting  rates  have  been  corrected  for  temperature  and,  more 
importantly,  for  RTG  background  (Appendix  I).  A  comparison  with 
the  Deep  River  Neutron  Monitor  data  (daily  means)  [Solar-Geo¬ 
physical  Data,  1572-1984]  is  shown  in  Figure  5.  A  plot  of 
monthly  mean  sunspot  numbers  [Solar-Geophysical  Data,  1984] 
is  given  in  Figure  6  to  provide  a  guide  to  solar  activity  for 
comparison  with  the  sample  of  Pioneer  10  data.  Some  general 
remarks  on  the  observed  data  are  as  follows: 

1.  Figures  2,  3,  and  4  show  that  solar  modulation  of  the 
cosmic  ray  intensity  extends  out  at  least  32  AU. 

2.  There  is  a  close  resemblance  in  considerable  detail 
between  the  intensity-time  curves  for  Pioneer  10  and  those  for 
Pioneer  11  throughout  the  past  eleven  years  (a  full 

solar  activity  cycle).  Identifiable  features  at  P10  are 
progressively  more  delayed  relative  to  the  corresponding 
features  at  Pll  as  the  difference  in  their  heliocentric 


distances  Ar  increases  [cf.  Van  Allen,  1 976].  The  delay  corre¬ 
sponds  roughly  to  the  ratio  of  Ar  to  440  km  s"1. 

3.  A  pronounced  26-day  periodic  variation  of  cosmic  ray 
intensity  with  amplitude  ~  4$  occurred  in  1974-76  (solar  activity 
minimum)  [Van  Allen,  1976]  but  such  periodicity  has  been  diffi¬ 
cult,  if  not  impossible,  to  discern  thereafter.  The  oblitera¬ 
tion  of  the  26-day  variation  is  attributed  to  an  increasing 
number  of  active  centers  on  the  sun  and  to  an  increasing  level 

of  roughness  in  the  interplanetary  magnetic  field. 

4.  The  correlation  between  the  interplanetary  intensity  of 
cosmic  rays  and  the  counting  rate  of  the  high  latitude  Deep 
River  Neutron  Monitor  (Figure  5)  is  impressively  good  (allowing 
for  propagation  delay)  despite  the  fact  that  the  effective 
primary  energy  to  which  the  latter  responds  is  several  times 
greater  than  the  effective  energy  to  which  our  GM  tubes  respond. 

A  detailed  comparison  for  two  great  Forbush  decreases  in  mid- 1982 
is  given  in  Figure  7. 

5>  In  Figure  6  there  is  seen  to  be  a  good  general  anti¬ 
correlation  between  cosmic  ray  intensity  and  sunspot  number, 
the  latter  being  taken  as  a  crude  measure  of  solar  activity. 
During  the  period  1972-1976  (declining  solar  activity)  the 
two  curves  are  approximately  in  phase  (±  1  month)  but  there¬ 
after  the  cosmic  ray  intensity  lags  behind  the  sunspot  number 
by  an  increasing  amount  (increasing  solar  activity)  and 


continues  to  do  so  during  the  subsequent  period  (maximum  and 
declining  solar  activity).  During  1980-1984  the  lag  is  about 
1.1  year  as  illustrated  by  the  portion  of  the  sunspot  curve  that 
has  been  moved  to  the  right  by  this  amount  and  superimposed  on 
the  cosmic  ray  intensity  curve  at  the  right  hand  end  of  the  latter 
in  Figure  6.  The  1.1  year  delay  is  much  greater  than  the  solar 
wind  propagation  delay  from  1  to  30  AU,  namely  ~  110  days. 
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III.  RADIAL  GRADIENT  OF  INTENSITY 


Making  the  common,  though  perhaps  incorrect,  assumption  of 
spherical  symmetry  about  the  sun,  the  omnidirectional  cosmic  ray 
intensity  J  is  a  function  of  two  variables  —  radial  distance  r  and 
time  t.  The  detailed  temporal  tracking  of  the  intensities  at  P10 
and  Pll  with  propagation  delay  At  «  (r^Q  -  r1^)/vsw  suggests 
the  further  assumption  that  the  radial  and  temporal  dependences 
of  J  are  uncoupled  and  that  J  can  be  written  as  the  product  of 
two  separate  functions,  viz. 


J10  (r10,  t  +  At)  =  K10  f(r10)  .  h(t  +  At)  (1) 

and 

Jn  (rir  =  K11  f(rn^  '  (2) 


wherein  K1Q  and  are  reciprocal  geometric  factors  (reciprocal 
effective  areas)  of  the  respective  detectors, 


At 


(r 


10 


rll} 


sw 


(3) 


and  vsw  is  the  mean  solar  wind  velocity  over  At. 
Then  by  (l)  and  (2) 
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J10  (r10’ t  +  At)  K10  f(r10)  h(t+At)  (4) 

Jn  (rii  ~  Kn  I^7uJ  * 


The  delayed  correlation  process  (see  next  section)  establishes 
At  and,  in  effect,  makes  equal  to  unity,  so  that 


(r 

11 


10 

Tr 


11 


t  +  At) 
>  t) 


K 


f<V 

^17 


(5) 


■where  K  s  %o^ll‘  ^  specific  assumption  on  the  form  of  f(r)  is 


f(r)  =  exp  (G  r)  . 


(6) 


Then 


10  (r10’ 
11  (rll 


t  +  At) 

TtT 


K  exp  (G  A  r)  . 


(7) 


G  is  called  the  integral  radial  gradient,  usually  expressed  in 
percent  per  AU.  The  scale  length  of  the  radial  dependence  is  l/G. 
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IV.  ANALYSIS  OF  DATA 

The  observed  data,  as  represented  by  Figures  2,  3,  and  4, 

have  been  analyzed  to  determine  the  dependence  of  intensity  on  r. 

First,  each  25-day  segment  of  daily  mean  counting  rates  of  a 

detector  on  P10,  say  C^Q,  is  compared  with  the  simultaneous  25-day 

segment  of  daily  mean  counting  rates  of  the  corresponding  detector 

on  Pll,  in  this  case  C^.  The  ratio  C^/C^  is  computed  for  each 

pair  of  days ;  then  the  25-day  mean  ratio  and  the  standard  deviation 

of  the  25  daily  mean  ratios  from  that  mean  are  computed.  The  process 

is  then  repeated  by  sliding  the  P10  block  of  days  backwards  in  time 

relative  to  the  Pll  data  in  increments  of  one  day  to  cover  the  range 

of  delays  corresponding  to  250  <  v  <  750  km  s-^.  The  case  that 

yields  the  lowest  value  of  the  standard  deviation  of  the  25  daily 

mean  ratios  from  the  25-day  mean  ratio  is  adopted  as  establishing 

the  propagation  delay  At  and  the  corresponding  25-day  mean  ratio. 

The  latter  quantity  is  f (r^)/f(r^^)  in  equation  (5)-  The  delays 

determined  by  the  above  procedure  are  not  of  as  high  an  accuracy  as 

was  possible  during  1974-76  [Van  Allen,  1976]  but  do  resemble 

Ar/v  with  a  typical  value  v  ~  440  ta  s  \  The  "delayed  ratios" 
sv  sv/ 

of  B^/B^,  ^'loA'll’  and  D]_c/Dll  are  Pitted  as  a  function  of  Ar  in 
Figures  8,  9,  and  10.  The  error  bars  on  the  ratios  are  standard 
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deviations  based  on  statistics  only,  with  no  contribution  from 
the  uncertainty  in  selecting  the  proper  pair  of  intervals  for 
which  the  ratio  is  computed. 

The  least  squares  linear  fits  to  these  semilog  plots  are 
of  the  form  exp  (GAr).  The  numerical  values  of  G  and  their 
formal  standard  errors  are  as  follows: 

For  +2.00  (±  0.07)  percent  per  AU 

For  C1Q/C11:  +2.07  (±  0.05) 

For  D10/D11:  +2.12  (±  0.0 6) 

The  three  overall  trend  lines  (Figs.  8,  9,  and  10)  for 
the  ratios  of  counting  rates  of  the  three  independent  pairs  of 
detectors  have  essentially  identical  slopes.  The  mean  value 
of  the  integral  radial  gradient  (Ep  >  80  MeV)  with  an  estimated 
uncertainty  from  all  known  causes  is 

G  =  +2.06  (±  0.20)  percent  per  AU,  1  <  r  <  32  AU.  (8) 

This  mean  radial  gradient  is  one  of  the  principal  results 
of  the  present  paper.  It  is  derived  from  data  over  an  eleven-year 
period.  The  calculation  of  the  gradient  is  based  on  the  varying 
quantity  Ar  whose  greatest  value  up  to  1  April  1984  is  16-7  AU. 

But  because  of  the  assumed  form 


f(r) 


exp  (Gr)  , 


(6) 
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the  gradient  G  has,  in  effect,  been  measured  piecewise  over  the 
range  1  <  r  <  32  AU. 

The  scale  size  of  the  heliospheric  modulation  region  is 
the  reciprocal  of  G,  namely 


50  AU. 


(9) 


V.  DISCUSSION 


Apart  from  two  regions  of  low  apparent  gradient  around  4r  w  2 
and  rj  8  AU  (Figs.  8,  9,  and  10),  the  value  of  G  has  been  fairly 
uniform  over  the  range  1  <  r  <  32  AU  and  over  a  solar  cycle  varia¬ 
tion  of  intensity  by  a  factor  of  2.7  (normalized  to  same  r).  This 
feature  of  similar  data  has  been  exploited  by  Webber  and  Lockwood 
[1981]  to  yield  an  estimate  of  the  radius  of  the  heliospheric  modula 
tion  region  ^  65  AU.  The  same  authors  using  data  for  >  60  MeV 
from  Pioneers  10  and  11,  Voyagers  1  and  2,  and  IMP'S  7  and  8  find  an 
overall  gradient  for  1  <  r  <  23  AU  of  +2.85  (±  0.5)  percent  per 
AU,  significantly  but  not  markedly  different  than  our  value  of 
+2.06  (±  0.020).  Also  Venkatesan  et  al.  [1984a]  using  independent 
data  from  Voyagers  1  and  2  (Ep  >  70  MeV)  and  IMP  8  (E  >  35  MeV) 
find  an  overall  gradient  for  1  <  r  <  13  AU  of  2  to  4  percent 
per  AU  and  favor  a  value  of  +3-5-  McKibben  et  al.  [1982]  using 
a  still  different  data  set  from  IMP's  5>  7,  and  8  and  Pioneers  10 
and  11  for  >  67  MeV  find  integral  radial  gradients  in  the  range 
+1.7  (±  0.6)  to  +3.9  (±  1.0)  percent  per  AU  for  various  periods  of 
time  and  various  combinations  of  data  from  different  spacecraft. 
These  authors  estimate  the  radius  of  the  heliosphere  as  50  to  70  AU 
at  both  solar  maximum  and  solar  minimum. 
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Honce  there  are  mild  disagreements  among  the  various  groups 
in  measuring  essentially  the  same  quantity  but  a  range  of  2  to  5 
percent  per  AU  for  the  integral  gradient  (E  >  60,  80  MeV)  now 
appears  to  be  reliable  and  is  perhaps  as  accurate  as  the  complex 
physical  situation  permits. 

McDonald  et  al.  [1981a]  using  data  from  Pioneer  10  and 
Helios  1  and  2,  the  latter  in  orbits  interior  to  that  of  the 
earth,  report  a  radial  gradient  out  to  22  AU  of  +3-5  percent  per 
AT J  in  the  differential  energy  range  115  <  <  220  MeV  and 

+3-6  (±  2)  to  +4.1  (±2)  percent  per  AU  for  helium  nuclei  in  the 
differential  specific  energy  range  114  <  eQ  <  580  MeV/nucleon. 

For  a  threshold  Ep  =  80  MeV  and  the  1965  (solar  minimum) 
proton  and  helium  spectra  of  Gloeckler  and  Jokipii  [1967]  our 
integral  gradient  corresponds  to  the  differential  gradient  at  an 
effective  energy  [Van  Allen,  1976]  in  the  range  420  to  890  MeV, 
depending  on  the  energy  dependence  of  the  gradient. 

Interpretation  of  a  radial  gradient  usually  follows  the 
heliospheric  line  of  thought  of  Morrison  [1956]  concerning  the 
11-year  modulation  and  Forbush  decreases  of  the  cosmic  ray  inten¬ 
sity.  In  early  work,  Parker  [1956]  advocated  a  geospheric  model 
for  both  of  these  effects  but  later  adopted  the  heliospheric 
model  and  led  the  way  in  giving  it  a  quantitative  structure 
[Parker,  1963]. 
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In  its  simplest  form  the  Morrison-Parker  model  contemplates 
the  emission  of  ionized  magnetized  gas  from  the  sun  with  originally 
no  cosmic  rays  therein.  As  the  gas  (plasma)  moves  outwards  through 
the  solar  system,  cosmic  rays  diffuse  into  the  gas  clouds  which  con¬ 
tain  "scattering"  centers  of  chaotic  magnetic  fields.  These  pre¬ 
sumed  scattering  centers  are  convected  outwards  from  the  sun, 
tending  to  carry  the  population  of  cosmic  rays  with  them.  In 
a  quasi-steady  state,  the  net  inward  diffusive  flux  of  cosmic  rays 
is  taken  to  be  equal  to  the  outward  convective  flux,  viz. 


n  v 


sw 


(10) 


where  n  is  the  number  density  of  cosmic  rays,  D  is  the  macro¬ 
scopic  radial  diffusion  coefficient,  and  vgw  is  the  solar  wind 
velocity.  Equation  (10)  ignores  factors  of  the  order  of  unity 
and  other  refinements  that,  though  possibly  significant,  decloud 
the  basic  issue  and  transcend  the  present  state  of  observational 
knowledge. 

It  is  further  supposed  (perhaps  incorrectly)  that  the 
physical  situation  is  spherically  symmetrical  and  that  the 
modulation  region  around  the  sun  is  spherical  with  outer  radius 
r  .  The  spherical  shell  bounding  this  region  is  called  the 

D 

modulation  boundary  (for  cosmic  rays).  Inasmuch  as  it  is 


impossible  for  the  sun's  influence  to  extend  infinitely  far,  a 
basic  question  is  the  value  of  r  .  In  the  spirit  of  (10),  the 
modulation  boundary  is  coincident  with  the  heliospheric  boundary, 
defined  as  the  spherical  shell  at  which  directed  flow  of  the 
solar  wind  ceases. 

Observationally,  passage  through  the  boundary  of  the 
heliosphere  (a  doubtless  fluctuating  quantity)  into  the  inter¬ 
stellar  medium  will  be  easily  recognized,  viz.: 

(a)  The  cosmic  ray  intensity  will  become  constant,  not  subject 
to  solar  modulation. 

(b)  The  radial  flow  of  the  solar  wind  will  cease. 

(c)  The  solar  wind  will  be  replaced  by  a  relatively  steady 
interstellar  wind,  having  a  much  different  velocity  vector, 
temperature,  composition,  etc. 

(d)  There  will  be  a  substantial  transition  in  the  vector 
magnetic  field. 

The  solar  system  is  in  motion  relative  to  a  selected  set 
of  stars  whose  proper  motions  and  radial  velocities  have  been 
measured.  From  these  measurements  it  is  inferred  that  the 
velocity  vector  of  the  solar  system  relative  to  the  "local 
standard  of  rest"  is  19-5  km  s  ^  toward  right  ascension  271°, 
declination  +50°  (the  "solar  apex”)  [Abell,  1982]  [Allen, 

1973].  In  ecliptic  coordinates  the  solar  apex  is  at  X  =  271°, 
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More  recently  and  of  greater  significance,  the  velocity  of 

the  solar  system  relative  to  the  immediately  nearby  interstellar 

medium  (consisting  in  part  of  neutral  atomic  hydrogen  and  helium) 

has  been  measured.  This  velocity  is  «  22  km  s  1  [Adams  and  Frisch, 

1977]  toward  right  ascension  252°,  declination  -17°,  or  toward 

\  =  2540,  3  =  +5°  [Ajello,  1979 ]•  The  medium  is  characterized  by 
1  -3  4 

n  «  0.04  -  0.08  cm  ,  T  «  10  K,  and  state  of  ionization 
H 

nH/nH+  »  0.5  [Pares ce,  1984]. 

These  facts  have  spawned  a  variety  of  speculation  that  the 
heliospheric  boundary  may  be  markedly  asymmetric.  But  because 
the  motion  is  at  a  velocity  of  only  22  km  s  1  whereas  the  solar1 
wind  velocity  is  »  440  km  s'1,  it  is  not  obvious  that  an  asymmetry 
is  significant,  especially  because  the  r.ra.s.  thermal  velocity 
(«  16  km  s'1)  of  the  atomic  hydrogen  and  the  sound  velocity  in  the 
interstellar  medium  are  comparable  to  the  relative  bulk  velocity. 

It  may,  nonetheless,  be  objected  that  our  determination  of 
G  by  comparing  data  from  P10  with  those  from  Pll  is  basically 
faulty  because  P10  is  moving  in  the  upstream  direction  relative 
to  the  interstellar  medium  whereas  Pll  is  moving  in  the  downstream 
direction  (Figure  l).  If  the  heliosphere  is  markedly  asymmetric 
along  that  line,  being  more  extended  downstream,  then  our  value  of 
G  would  be  falsely  low.  Other  work,  however,  comparing  both  P10 
and  Pll  data  out  to  a  common  distance  of  12  AU  with  IMP  data  near 
the  earth  [Webber  and  Lockwood,  1981]  gives  no  evidence  for  such 


asymmetry. 


The  time  dependences  of  the  heliographic  latitudes  of  P10 
and  Pll  are  shown  in  Figure  11.  In  previous  work  [Van  Allen, 
1576]  on  the  data  for  1975*  we  found  a  latitudinal  gradient  of 
0  (±  2)  percent  per  10°  of  heliographic  latitude.  Other  recent 
attempts  to  determine  the  latitudinal  gradient  have  been  made  by 
Roelof  et  al.  [1983],  Decker  et  al.  [1984],  and  McKibben  et  al. 
[1979]-  No  consensus  has  yet  emerged  nor  is  one  likely  before 
high  latitude  observations  by  the  International  Solar  Polar 
Mission  in  1986  ff. 

For  many  years,  Lockwood  [i960]  has  urged  the  view  that 
the  11-year  solar  cycle  of  cosmic  ray  intensity  is  caused  by  a 
higher  incidence  of  episodic  depressions  of  the  intensity 
(Forbush  decreases)  during  solar  maximum  activity  from  which 
recovery  is  sufficiently  slow  that  such  overlapping  decreases 
have  a  cumulative  effect.  A  convincing  case  for  this  view  is 
given  in  Lockwood  and  Webber  [1984].  We  are  in  agreement  with 
this  explication  as  is  readily  apparent  from  Figures  2,  3,  and  4, 
and  earlier  work  [Van  Allen,  1979]-  In  Figure  7  we  give  a 
detailed  presentation  of  P10,  Pll,  and  Deep  River  Neutron 
Monitor  data  for  two  noteworthy  Forbush  decreases  during  1982. 
Provisionally,  we  identify  the  effects  on  DOY's  l6l  at  the 
earth,  187  at  Pll,  and  215  at  P10  as  attributable  to  the  . jme 
change  in  interplanetary  conditions  and  the  effects  on 
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DOY's  195,  219,  and  271  as  attributable  to  another.  The  two 
episodes  are  different  in  detail  but  both  resulted  in  a  durable 
depression  of  intensity  at  both  P10  and  Pll  [cf.  Van  Allen, 
1979].  The  details  of  the  two  cases  are  given  in  Table  1. 

From  the  definition  of  G  and  equation  (10), 


1  Sn  sw 

G  "  n  5F  “  ~  * 


(11) 


-1 


Using  observed  values  of  G,  0.02  (AU)  ,  and  v  , 


440  bn  s  \ 


D  =  3-3  X  1022  cm2  s-1  =  4600  (AU)2  yr_1.  (IP) 


The  corresponding  scattering  mean-free-path, 

=  =  3.5  x  1012  cm  =  o.22  AU.  (13) 

c 

This  value  of  the  mean-free-path  may  be  compared  to  the  gyro- 
radius  of  a  proton  of,  say,  2  GeV  energy  in  an  interplanetary 
magnetic  field  of  X  10~^  gauss,  namely  0.06  AU. 

The  above  calculations  are  intended  only  to  fix  the  orders 
of  magnitude  of  quantities  and  make  no  pretense  to  accuracy  (e.g. , 
D  may  be  a  function  of  r,  t,  and  particle  energy). 


!*&<• 
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We  now  have  a  direct  measurement  that  r  >  32  AU  and  a 

B 

reliable  estimate  of  the  scale  length  of  the  modulation  region, 
namely  50  AU.  Nonetheless,  r  may  well  be  considerably  greater 
than  50  AU. 

One  simple  method  of  estimating  r  is  to  suppose  that  the 

B 

boundary  of  the  heliosphere  is  at  the  radial  distance  at  which 
the  dynamic  pressure  of  the  solar  wind  equals  the  interstellar 
magnetic  field  pressure  B  /8rr,  namely, 

nmv2  »  B2/8tt  .  (l4) 

sw  '  '  ' 

Taking  n  =  5/r  ,  vgw  =  440  km  s"  ,  m  =  1.67  X  10"  g,  and 

B  =  7  x  10  6  gauss  [Allen,  1973],  one  finds  from  (14)  that 

rB  »  90  AU.  But  the  uncertainties  in  the  various  quantities 

do  not  permit  much  confidence  in  the  result  nor  do  they  permit 

much  confidence  in  claiming  a  variation  of  r  with  solar 

B 

activity. 

Another  method  of  estimating  r^  is  to  use  modulation 
theory,  such  as  it  is,  to  "demodulate"  the  spectrum  to  inter¬ 
stellar  space  and  then  to  estimate,  using  observed  values  of 
the  differential  gradient  g,  the  value  of  r  at  which  the 
demodulated  spectrum  would  occur.  By  this  means  McKibben  et  al. 

[1982]  estimate  that  r  is  50  to  70  AU,  at  both  solar  maximum 

B 


and  solar  minimum. 
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We  have  noted  above  that  Webber  and  Lockwood  [1961]  have 

used  a  kindred  method  to  estimate  that  r  >  65  AU. 

B 

Still  another  method  of  estimating  r  is  to  suppose  that  the 

B 

volume  density  of  cosmic  ray  kinetic  energy  in  interstellary  space 
is  equal  to  B  /8tt  therein.  Such  an  "equipartition"  of  energy  is 
discussed  critically  by  Cowling  [1967],  who  gives  some  support  to 
that  concept  but  also  expresses  doubt  as  to  its  validity.  None¬ 
theless  it  yields  results  of  some  interest.  Using  the  proton  and 
helium  spectra  of  Gloeckler  and  Jokipii  [1967],  we  calculate  the 
■volume  density  6  of  galactic  cosmic  rays  at  1  AU  at  solar  minimum 
by  numerical  integration  of  the  following 


l 


E 

3 


(15) 


wherein  3  is  the  ratio  of  the  speed  of  the  particle  to  the  sp.-od 
of  light  c  and  dj/dE  is  the  differential  unidirectional  spectral 
intensity.  The  results  are 


e 

p 


7.3  X  10  15  erg  cm  ^ 


(about  16$  of  which  is  for  particle  energies  below  ?()00  MeV,  tlv- 
modulation  sensitive  energy  region) 


ami 


&a  =  3.2  x  10  ^  erg  cm  ^ 


2b 


(about  19$  of  which  is  below  2000  MeV/nucleon) .  The  sum 

e  =  £  +  £  =  1.0  x  10"12  erg  cm"5.  (l6) 

pa  '  ' 

-6 

For  an  interstellar  magnetic  field  B  =  7  X  10  gauss  [Allen 

1973] 

2 

=  1.9  X  10-12  erg  cm"5  .  (17) 

Thus  it  appears  possible  for  the  modulation-sensitive  portion  of 

2 

the  spectrum  to  increase  by  a  factor  of  «  e  above  its  1  AU  value 
in  order  to  approximate  qua si-equi partition  in  the  interstellar 
medium.  This  line  of  thought  suggests  that  r0  «  (2)(50)  =  100  AU. 
Finally,  it  appears  possible  to  base  an  estimate  of  r^  on 

1J 

the  1  year  lag  in  response  of  the  cosmic  ray  intensity,  even  at 
«  30  AU,  to  solar  activity  as  measured  by  sunspot  numbers 
(Figure  6)  [cf.  FuLks,  1575l>  A  grossly  simplified  model  for 
this  effect  is  along  the  following  lines.  First  suppose  that  a 
steady  state  has  been  established  with  a  particular  value  of  D 
throughout  the  heliosphere.  Then  suppose  that  the  conditions  at 
the  sun  change  suddenly  so  that  D  increases  by,  say,  a  factor  of  2. 
The  spherical  region  of  higher  D  expands  outward  at  velocity 
v  until  it  reaches  the  heliospheric  boundary,  wherever  that 
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may  be.  The  time  constant  t  for  the  increase  of  intensity  interior 
to  the  expanding  region  is  of  the  order  of 


T 


«  0. 5 


(18) 


Hence 


r 


B 


D  t 

0.5 


2  -1 

’wherein  D  «  4600  (AU)  yr  and  t  »  1  yr.  By  (19), 


(19) 


r  «  100  AU  . 

B 

In  one  year,  the  solar  wind  travels  «  93  AU,  giving  some  plausi¬ 
bility  to  this  result. 

None  of  the  foregoing  estimates  of  r„  withstands  critical 
scrutiny  but  they  do  give  some  feeling  for  the  nature  of  the 
problem. 

Observationally,  r,,  >  32  AU  and  it  would  not  be  astonish- 
B 

ing  if  its  value  lies  between  50  and  100  AU.  The  value  of  r  may 

B 

systematically  decrease  during  solar  minimum  conditions  and  in¬ 
crease  during  solar  maximum  conditions  and  may  moreover  undergo 
sporadic  fluctuations  of  many  AU  as  the  dynamic  pressure  of  the 
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solar  -wind  fluctuates.  Observations  during  the  next  decade 
should  yield  an  improved  estimate  of  rB  or  perhaps  an  actual 
determination. 
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APPENDIX  I.  Remarks  on  Validity 
and  Accuracy  of  Data 


1.  Individual  counting  rates  of  a  number  of  miniature  Geiger- 
Muelier  tubes  provide  the  basic  data  for  the  present  paper.  These 
GM  tubes  were  built  by  the  EON  Corporation.  All  of  them  are  halogen- 
quenched,  "infinite  lifetime"  tubes.  Detectors  A,  B,  C,  and  G  on 
Pioneer  10  (Van  Allen  et  al.,  1974)  and  A,  B,  and  C  on  Pioneer  11 
(Van  Allen  et  al. ,  I960)  are  EON  Type  6213.  Similar  tubes 
have  been  used  previously  in  University  of  Iowa  instruments  on 
Explorer  4,  Pioneer  3,  Pioneer  4,  Explorer  7,  Injun  1,  Explorer  12, 
TRAAC,  Mariner  2,  Explorer  14,  Explorer  15,  Injun  3,  0G0  1, 

Injun  4,  Mariner  4,  000  2,  Explorer  33,  Mariner  5,  Explorer  35, 

Injun  5,  and  Hawkeye  1.  Tubes  in  those  instruments  logged  over 
700,000  detector-hours  of  reliable  operation  in  flight.  Many 
individual  tubes  experienced  over  10^  counts  and  a  few,  over 
10'1''*"  counts. 

Detectors  D,  E,  and  F  on  each  spacecraft  are  BON  Type  5107, 
being  of  different  size  and  shape  than  the  6213 's  as  described  in 
detail  in  the  second  of  the  foregoing  references.  Tubes  of  this 
type  had  not  been  flown  before  but  were  found,  after  selection 
and  testing,  to  have  excellent  reliability. 


The  tubes  that  were  assembled  in  the  flight  units  were  the 
survivors  of  a  rigorous  selection  procedure  on  a  large  batch  of 
each  of  the  two  types  of  tube.  The  procedure  involved  mechanical 
vibration  testing  and  then  running  the  plateau  curves  (counting 
rate  vs.  applied  voltage  when  exposed  to  a  fixed  7-ray  source) 
as  a  function  of  temperature  and  of  exposure  age.  Of  the  order 

Q 

of  10°  counts  were  accumulated  on  each  tube.  About  10 %  of  the 
tubes  were  judged  suitable  for  flight.  The  selected  tubes  were 

O 

assembled  in  flight  units  and  again  subjected  to  ~  10°  counts  by 

60 

hard  x-ray  and  Co  sources  for  the  determination  of  apparent  vs. 
true  counting  rate  curves  before  and  after  an  extended  program  of 
thermal -vacuum  testing.  Each  of  the  completed  instruments  includ¬ 
ing  detectors  operated  properly  over  the  temperature  range  -20° 
to  +100°  F.  The  temperature  coefficient  of  counting  rate  for  a 
fixed  Co^°  7-ray  source  was  determined  for  each  tube  (and  its 
associated  circuitry)  as  assembled  in  the  flight  unit  over  the 

temperature  range  -20°  to  +90°  F.  For  Pioneer  10,  the  mean 

-4 

coefficient  was  3-98  X  10  per  degree  Fahrenheit  for  A,  B,  C, 

-4 

and  G  and  -2.18  X  10  per  degree  Fahrenheit  for  D,  E,  and  F. 

.  -4 

For  Pioneer  11,  the  mean  coefficient  was  3.69  x  10  per  degree 
Fahrenheit  for  A,  B,  and  C  and  -6.5  X  10~^  for  D,  E,  and  F. 

All  observed  flight  rates  are  corrected  to  +75°  F  using 
the  above  coefficients  and  the  measured  temperature  of  the 
detector  assembly  (see  next  paragraph). 
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The  in-flight,  temp.  ■railin'  of  a  calibrated  thermistor 
.imbedded  in  our  detector  assembly  is  telemetered  continuously. 

For  both  Pioneer  10  and  Pioneer  11,  the  highest  temperature 
(+77°  F)  occurred  on  the  day  of  launch.  The  temperature  declined 
rapidly  to  about  48°  F  during  the  first  two  months  and  then  more 
slowly  to  +40°  F  (P  10)  and  +45°  F  (P  11)  b^  the  end  of  the  first 
year  of  flight.  Thereafter  the  temperatures  have  declined  linearly 
with  time  at  the  rate  of  2 0  F  per  year.  As  of  April  1984,  both 
temperatures  were  +20.6°  F  (following  turn-off  of  power  to  the  IPP 
instrument  on  Pioneer  11).  These  temperatures,  essentially  inde¬ 
pendent  of  the  intensity  of  sunlight  for  the  past  ten  years,  are 
dependent  principally  on  the  power  dissipated  in  the  instrument 
(0.72  watt)  and  in  other  instruments  on  the  instrument  platform 
of  the  spacecraft.  We  anticipate  no  problem  with  temperature  during 
the  prospective  lifetime  of  either  spacecraft.  Power  sharing 
(i.e.,  alternate  on-off  operation)  among  the  instruments  on 
Pioneer  11  will  be  initiated  in  late  1984,  and  on  Pioneer  10 
about  two  years  later.  This  procedure  will  lower  the  "power-off" 
temperature  of  our  detector  assembly  by  ^  8°  F. 

Because  of  the  gradual  and  small  change  of  the  temperatures 
during  most  of  the  duration  of  the  two  missions  and  of  the  small 
magnitude  of  the  temperature  coefficients,  the  estimated  uncer¬ 
tainties  in  the  adopted  values  of  the  coefficients  have  a  trivial 
effect  on  the  results  of  the  present  study. 


30 


3.  The  spacecraft  bus  voltage  is  regulated  at  28.0  V  (±  V%) 
and  the  terminal  voltage  of  the  internal  power  supply  for  the 
electronics  in  our  instrument  is  separately  regulated  at 
+7-92  V  (±  2 i).  Both  of  these  voltages  are  telemetered  con¬ 
tinuously  via  an  ana log-to -digital  converter  having  64  linear 
steps.  Both  voltages  on  both  spacecraft  have  been  invariant 
throughout  their  flight  periods.  A  change  of  one  digital  step 
would  be  a  matter  of  some  concern  but  our  instrument  is  unaffected 
by  a  change  in  bus  voltage  over  the  range  26  to  29  V. 

4.  No  one  of  the  GM  tubes  has  exhibited  any  anomaly  (i.e., 
spuriously  high  or  low  counting  rates  or  failure  to  count)  at 
any  time  before  or  since  launch  (now  over  twelve  years  of  flight 
for  Pioneer  10,  including  passage  through  the  outer  and  inner 
radiation  belts  of  Jupiter,  and  over  eleven  years  of  flight  for 
Pioneer  11,  including  passage  through  the  outer  and  inner  radia¬ 
tion  belts  of  both  Jupiter  and  Saturn).  Each  of  the  least 
shielded  GM  tubes  (Detector  A)  has  experienced  about  3  X  10^° 
counts  during  the  missions  to  date  and  the  more  heavily  shielded 
tubes,  about  an  order  of  magnitude  fewer  counts.  The  dominant 
contribution  occurred  during  passage  through  Jupiter's  magneto¬ 
sphere  in  both  cases.  No  discontinuous  change  in  interplanetary 
counting  rates  of  any  detector  before  and  after  the  encounters 


was  observed. 


Inasmuch  as  it  was  not  feasible  to  provide  an  in-flight 
radiation  source  for  an  end-to-end  check  of  the  GM  tubes  and 
their  associated  electronics,  there  is  no  direct  basis  for  ruling 
out  gradual  drifts  in  sensitivity.  There  are,  however,  four 
grounds  for  considering  drifts  unlikely: 

(a)  Extensive  laboratory  experience; 

(b)  The  likelihood  that  any  such  drift  would  be  character¬ 
istic  of  an  individual  tube. 

(c)  The  accurate  constancy  of  the  ratios  of  the  interplanetary 

counting  rates  among  the  several  tubes  on  the  same  space¬ 
craft,  e.g. ,  etc.  throughout  the  missions;  and 

(d)  The  close  agreement  of  the  time-variability  of  the  three 

independent  ratios  of  the  counting  rates  of  correspond¬ 
ing  tubes  on  the  separate  spacecraft,  b;lo/b;q>  ^lO^ll' 
and  (see  text). 

Points  (c)  and  (d)  support  (b)  by  the  additional  fact  that,  as  noted 
above.  Detectors  D  are  of  a  quite  different  type  than  are 
Detectors  B  and  C. 

5«  Detector  A  on  Pioneer  11  has  no  shielding  over  its  thin 
mica-window  and  is  sensitive  to  low  energy  particles;  it  detects 
numerous  solar  and  iriterplanetary  particle  events  so  that  its 
rate  is  seldom  at  the  cosmic  ray  background  level.  Also, 

Detector  G  on  Pioneer  10  is  a  GM  tube  but  such  a  tube  was 
replaced  by  a  solid  state  detector  on  Pioneer  11.  Hence, 
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in  the  present  study  we  have  settled  on  study  of  the  three 

ratios  B^/B^,  Cic/Cil’  and  Dic/Dll‘  11116  corresPonding1y 
designated  detectors  on  the  two  spacecraft  have  virtually 
identical  shielding  and  physical  placement  in  the  instrument 
package.  However,  we  find  that  virtually  the  same  results 
are  obtained  using  ratios  such  as  C^q/B^,  etc. 

6.  The  concordant  results  (see  text)  among  the  three  adopted 
counting  rate  ratios  provide  overall  confidence  in  the  validity 
of  the  observational  data,  irrespective  of  detailed  considera¬ 
tions  . 

7.  On  the  basis  of  theoretical  failure  analysis,  two  inde¬ 
pendent  signal  processors  (for  converting  the  GM  pulses  to  a 
quasi-logarithmically-scaled  digital  output)  were  included  in 
each  instrument.  During  Pioneer  11 's  encounter  with  Jupiter  in 
December  1974,  a  failure  (presumably  attributable  to  cumulative 
radiation  damage)  occurred  in  high  order  bits  of  one  of  the  two 
channels  of  the  "standby"  processor.  By  radio  command,  this 
processor  was  replaced  by  the  "main"  processor.  The  latter  has 
functioned  properly  throughout  the  subsequent  flight.  No  anomaly 
of  any  kind  at  any  time  has  occurred  in  either  of  the  two 
processors  on  Pioneer  10. 

8.  The  electrical  power  for  each  of  the  two  spacecraft  is 
provided  by  four  Type  SNAP  19  Radioisotope  Thermal  Generators 
(RTG's)  of  the  Atomic  Energy  Commission.  The  heating  of  the 
thermoelectric  system  is  provided  by  a  mixture  of  plutonium 


(half  life  8 6.4 


isotopes,  the  most  important  of  which  is  Pu 

258 

years).  The  total  activity  of  Pu  in  the  four  units  is 

77,000  curies  (2500  watts  thermal).  There  are  several  much 

weaker  Radioisotope  Heater  Units  (RHU’s),  also  using  plutonium, 

for  spot  heating.  The  typical  isotopic  composition  by  weight 

[Kaminskas  et  al.,  19711  of  the  plutonium  "fuel"  by  nuclear  mass 

number  is  238,  80$;  239,  15-8$;  240,  2.8$;  24l,  1$;  242,  0.4$; 

and  236,  1  part  per  million.  Basic  data  on  half-lives,  decay 

schemes,  etc.  are  known  for  all  of  these  isotopes.  Using  these 

data,  extensive  overall  data  on  the  RTG's,  and  laboratory  tests 

236 

on  our  apparatus,  we  have  concluded  that  Pu  ,  despite  its  low 

relative  abundance,  is  the  dominant  constituent  in  causing 

increasing  background  counting  rates  of  our  CM  detectors  as  a 

function  of  time  following  fuel  separation.  The  next  most 

important  contribution  (declining  with  time)  to  the  background 

234 

rates  comes  from  L  x-rays  of  energy  ~  100  keV  from  U  ,  the 

238 

long-lived  daughter  of  the  alpha  particle  decay  of  Pu 

The  complex  decay  chain  of  Pu2^  (half-life  2.85  years) 
and  the  emissions  from  successive  progeny  have  been  described 
and  discussed  by  Kaminskas  et  al.  [1971],  Thomsen  [1974],  and 
Thomsen  and  Van  Allen  [1976]  and  have  been  reviewed  recently 
by  the  current  authors.  The  principal  conclusion  is  that  the 
abundance  of  T A208  in  the  Pu256  decay  chain  is,  following  a  year 
or  two  of  aging,  the  most  important  measure  of  counting  rate 


background  by  virtue  of  the  2.6l  MeV  7-ray  from  an  excited  state 

POo  POft 

of  its  stable  daughter  Pb  .  The  abundance  of  T l  is,  in 

turn,  proportional,  at  any  given  time,  to  the  abundance  of  one  of 

ppQ 

its  ancestors  Th  .  The  latter  increases  from  essentially  zero 

at  the  time  of  fuel  separation  to  a  maximum  value  18.6  years 

later  and  then  declines  gradually. 

Prelaunch  tests  with  the  flight  RTG's  and  our  flight 

instrument  on  the  spacecraft  were  made  in  the  assembly  room  of 

TRW  Systems,  Inc.  for  the  purpose  of  determining  the  background 

counting  rates  of  all  GM  tubes.  Also,  simulated  aging  tests 
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were  done  by  adding  Th  sources  of  calculated  activity  during 
these  tests.  In  our  earlier  work  [Van  Allen,  1972;  Thomsen,  197*+; 
Thomsen  and  Van  Allen,  1976;  Van  Allen,  1979>  1980]  on  the  cosmic 

ppQ 

ray  gradient,  a  theoretical  curve  of  the  activity  of  Th  as  a 
function  of  time  was  fit  to  these  test  data  and,  despite  some 
mild  internal  inconsistencies,  this  curve  was  used  to  calculate 
the  subtractive  corrections  to  all  observed  counting  rates. 

9.  During  the  Saturn  encounter  of  Pioneer  11  there  occurred 
a  unique  and  precious  opportunity  to  check  the  magnitude  of  the 
RTG  corrections  under  free  space  conditions.  The  spacecraft 
passed  under  the  dense  Rings  A  and  B  of  Saturn  into  periapsis 
at  1.2*+  Rg.  The  mean  counting  rates  of  GM  tubes  A,  B,  and  C 
are  plotted  as  a  function  of  radial  distance  in  Figure  12. 

After  correction  for  the  contribution  of  primary  cosmic  rays 
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(using  St^rmer  cut-off  theory)  and  of  secondaries  produced  in 
the  rings  after  the  fashion  of  Chenette  et  al.  [1980],  we  find 
the  residual  mean  counting  rate  of  A,  B,  and  C  to  be 
R11  =  0*272  ±  0.010  c  s'1  (Figure  12).  This  rate  is  assigned 
solely  to  the  RTG's. 

Using  our  previous  formulae  based  on  preflight  tests 
a  rate  of  0.50  c  s  ^  was  calculated  for  this  time.  This  important 
discrepancy  is  not  really  understood  but  is  thought  attributable 
to  room  scattering  during  the  preflight  tests  and  uncertainties 
in  the  simulated  aging  tests.  We  now  believe  that  the  pre¬ 
viously  published  strong  upturn  in  the  radial  dependence  of 
Ciq/Ci^,  etc.  following  the  major  Forbush  decrease  of  cosmic  ray 
intensity  in  April-May  1978  is  spurious,  having  been  caused  by 
the  excessively  large  RTG  corrections  as  applied  to  the  then  much 
reduced  rates.  [See  Figure  4  (a),  (b),  and  (c)  of  Van  Allen, 

1980. ]  Previous  to  that  time  the  ratios  were  much  less 
affected  by  the  excessive  correction. 

We  have  adopted  the  Saturn  encounter  value  of  R^  and  be¬ 
lieve  it  to  be  reliable  to  ±  5$*  The  corrections  for  individual 
detectors  A,  B,  C,  and  D  are  derived  from  the  mean  value  in 
accordance  with  the  ratios  of  their  interplanetary  counting  rates 
during  quiet  conditions,  viz.: 
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|  R  (Au)  =  1.137  Rn 

R  (Bu)  =  0.905  Ru 
R  (Cu)  =  0.958  Rlx 
|  R  (Du)  =  0.254  Rlx 

Our  working  formula  for  calculating  the  time  dependence  of  the 

(  mean  subtractive  correction  R^  is  based  on  AEC  supplied  informa  - 

I 

tion  on  the  isotopic  composition  of  the  several  RTG’s,  dates  of 
separation  of  the  fuel,  etc.,  and  on  basic  radioactive  decay-chain 
data  (cf.  Thomsen  and  Van  Allen,  1976).  The  correction  curve  was 

* 

fit  to  the  Saturn-encounter  point  and  to  the  preflight  point  (as 
adjusted  downward  for  estimated  room  scattering).  The  curve  is 
.  principally  dependent  on  the  former,  well-determined  point  and 

only  weakly  dependent  on  the  re-estimated  preflight  point.  Our 
adopted  algebraic  formula  is 

i  -  2 

Ru  =  1.103  X  10  *  [2  f(t)  +  0.9  f  (t  +  0.832)] 

+  5.53  x  10”^  [3  g(t)  +  g  (t  +  0.832)]  c  s”1 

where 

‘  f  (t)  =  23.792  exp  (-  0.3623  t) 

+  12.120  exp  (-  0.0094  t) 

-  35-912  exp  (-  0.2432  t), 

|  g  (t)  =  exp  (-  0.0081  t), 

I 
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nnd  1.  in  measured  in  years  with  t  =  0  on  15  March  1971  (HOY  74). 

A  plot  of  R^  vs.  time  is  given  in  Figure  13- 

There  is  no  corresponding  determination  of  the  in-flight, 
free-space  value  of  the  RTG  corrections  for  Pioneer  10.  Hence 
we  have  assumed  that  they  can  be  derived  from  the  corrections  for 
Pioneer  11,  using  the  somewhat  different  values  of  isotropic 
composition  of  the  plutonium,  dates  of  separation,  etc.  Our 
adopted  algebraic  formula  is 

R10  =  3.862  X  10"2  f(t)  +  2.211  X  10"2  g(t)  c  s-1 

where  f(t)  and  g(t)  are  as  previously  specified  and  t  is  the  time 
in  years  with  t  =  0  on  15  June  1970  (DOY  166).  A  plot  of  R^q 
vs.  time  is  also  given  in  Figure  13. 
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10.  The  duty  cycle  for  accumulating  counts  from  each  of  the 
GM  tubes  A,  B,  C,  and  D  is  l/ll,  irrespective  of  telemetry 
bit  rate.  Thus,  in  a  solid  24  hours  of  telemetry  reception  at 
a  typical  counting  rate  (for  A,  B,  or  C)  of  0.8  c  s-1,  6300  counts 
are  recorded.  The  statistical  standard  deviation  of  such  a  one-day 
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mean  is  1.3$  or  of  a  25-day  mean,  0.3 $.  The  actual  telemetry 
reception  approached  100$  early  in  the  missions  but  later  has 
been  more  typically  in  the  range  8  to  16  hours  per  day  with  a 
corresponding  increase  in  relative  statistical  standard  devia¬ 
tion.  The  typical  counting  rate  of  the  GM  tubes  D  is  about 
0.3  c  s  and  the  statistical  standard  deviation  is  larger  than 
those  computed  above  by  a  factor  of  1.6. 

The  r.m.s.  fluctuation  of  mean  daily  rates  from  the  mean 
rate  for  a  25-day  period  is  never  less  than  several  times  as 
great  as  the  statistical  expectation -- presumably  signifying 
real  physical  fluctuations. 
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Table  1 

Data  on  Two  Large  Forbush  Decreases 
(Fig.  7) 


1982 

i 

r 

\ 

First  Event 

BOY 

Decrease 

AU 

Deep  River  Neutron  Monitor 

161  ±  1 

10 

1.02 

258  .°8 

Pioneer  11 

187  ±  2 

18 

12.1 

220.3 

Pioneer  10 

215  ±  1 

Ik 

27.9 

65.5 

Second  Event 

Deep  River  Neutron  Monitor 

195  ±  1 

21 

1.02 

291.2 

Pioneer  11 

219  ±  1 

21 

12.3 

221.5 

Pioneer  10 

271  ±  1 

16 

28.3 

65.7 

Notes : 

(1)  The  two  events  are  qualitatively  different  in  nature. 

(2)  The  apparent  velocities  of  propagation  Ar/At  (ignoring 
differences  in  longitude)  for  the  first  event  are 

7^0  ±  60  km  s  from  the  earth  to  Pll  and  860  ±  30  km  s"'1' 
from  the  earth  to  P10. 

(?)  The  apparent  velocities  of  propagation  for  the  second 
event  are  810  ±  50  km  s  ^  from  the  earth  to  Pll  and 
620  ±  15  km  s  ■*"  from  the  earth  to  P10.  The  discrepancy  in 
apparent  velocities  in  this  case  is  presumed  to  be  due  to 
the  more  directional  nature  of  the  second  event,  as  evidenced 
by  the  neutron  monitor  curve,  and  the  neglect  of  corotation 


time  to  reach  P10. 
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CAPTIONS  FOR  FIGURES 

Fig.  1.  Ecliptic  plane  projection  of  the  trajectories  of 

Pioneer  10,  Pioneer  11,  Voyager  1,  and  Voyager  2. 

Also  shown  are  the  orbit  of  the  earth  and  portions  of 
the  orbits  of  Jupiter,  Saturn,  and  Uranus;  and  a  table 
of  dates  of  launch  and  closest  approach  of  the  four 
spacecraft  to  Jupiter,  Saturn,  and  Uranus.  The  X-axis 
points  to  the  earth's  vernal  equinox  (\  =  0°  at  this 
point  and  increases  counterclockwise). 

Fig.  2.  Fully  corrected  daily  mean  counting  rates  of  GM 

detectors  B  on  Pioneers  10  and  11  as  a  function  of  time. 
Heliocentric  radial  distances  are  given  by  the  upper 
scale  for  Pioneer  11  and  the  lower  scale  for  Pioneer  10. 

Fig.  3.  Same  as  Figure  2,  except  for  GM  detectors  C. 

Fig.  4.  Same  as  Figure  2,  except  for  GM  detectors  D. 

Fig.  5-  A  repeat  of  Fig.  3  with  the  addition  of  the 

daily  mean  counting  rate  (of  the  Deep  River  Neutron, 

Lat  46.°10  N,  Long  282.°50  E,  alt  145  m).  The  ordinate 
is  in  counts  hr"V300  [Solar-Geophysical  Data,  15r7’2-1984]. 
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Fig.  6.  Comparison  plots  of  the  counting  rates  of  GM 

detector  C  on  Pioneer  10  and  monthly  mean  sunspot  numbers 
[Solar-Geophysical  Data,  1972-7*0-  Note  that  the  latter 
plot  is  inverted.  A  portion  at  the  right  hand  end  of 
the  sunspot  number  plot  has  been  moved  to  the  right 
by  1.1  year  and  superimposed  on  the  cosmic  ray  intensity 
plot  in  the  upper  panel  to  illustrate  the  time  delay. 

Fig.  7*  Graphical  comparison  of  the  data  from  Pioneer  10 

and  Pioneer  11  with  those  from  the  Deep  River  Neutron 
Monitor  for  two  large  Forbush  decreases  in  mid-1982. 

Fig.  8.  Semilogarithmic  plot  of  the  "delayed"  ratios  of 

twenty-five -day  mean  counting  rates  of  GM  detectors  B 
on  Pioneer  10  and  on  Pioneer  11  vs.  A r,  the  difference  in 
the  heliocentric  radii.  The  straight  line  is  the  least- 
squares  fit  to  the  data  of  an  equation  of  the  form 
exp  (GAr)  . 

Fig.  9>  Same  as  Figure  8,  except  for  GM  detectors  C. 

Fig.  10.  Same  as  Figure  8,  except  for  GM  detectors  D. 

Fig.  11.  Heliographic  latitudes  of  Pioneer  10  and  Pioneer  11 


vs.  time. 


Fig.  I..'1. 
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Moan  counting  rates  of  GM  detectors  A,  B,  and  C 
on  Pioneer  11  vs.  radial  distance  from  the  center  of 
Saturn  during  the  1  September  1979  encounter. 

Fig.  13-  Adopted  curves  of  the  time  dependence  of  the  mean 

RTG  corrections  to  the  counting  rates  of  04  detectors 
A,  B,  and  C  on  Pioneer  10  and  Pioneer  11.  The  corrections 
are  subtracted  from  the  raw  counting  rates  after  tempera¬ 
ture  correction  of  the  latter. 
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